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Abstract 

In this paper we deal with the dimension of multisequences and related properties. 
For a given multisequence W and R € Z+, we define the R— extension of W. Further we 
count the number of multisequences W whose R— extensions have maximum dimension 
and give an algorithm to derive such multisequences. We then go on to use this theory 
to count the number of Linear Feedback Shift Register (LFSR) configurations with multi 
input multi output delay blocks for any given primitive characteristic polynomial and 
also to design such LFSRs. Further, we use the result on multisequences to count the 
number of Hankel matrices of any given dimension. 

1 Introduction 

Linear recurring sequences, over a finite field ¥ q , with maximum period have been shown 
to exhibit several important randomness properties such as 2-level autocorrelation property 
and span-n property (all nonzero subsequences of length n occur once in every period) [T]. 
As a result they find applications in a wide array of areas including cryptography [2], error 
correcting codes [3] and spread spectrum communication [I]. 

An obvious extension of a sequence of scalars is a sequence of vectors over the given finite 
field. Such a sequence of vectors is known as a multisequence. Periodic multisequences and 
linear relations among the elements of such sequences, have been a subject of study for a 
considerable period of time [5], [6], [7]. Generating multisequences with a given minimal 
polynomial has been an important problem motivating papers like [5], [H] and |10j . 

In this paper, we start by deriving some basic theorems regarding multisequences. We 
then introduce the concept of an R— extension of a multisequence, for R G Z™. We then 
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derive an algorithm to generate multisequences whose R— extensions have maximum dimen- 
sion. Further we derive a formula for the number of multisequences having this property. As 
an application, we show that the problem of generating some special class of LFSR configura- 
tions for any given primitive characteristic polynomial is a special case of the above problem. 
We then go on to count the number of such LFSR configurations using the formula derived 
for multisequences. Finally we demonstrate another application of the theory developed : a 
novel way to count the number of full rank Hankel matrices with entries from a given finite 
field. 

In the remainder of this paper, ¥ q denotes a field of cardinality q, where q is a prime 
power. ¥ q [s] denotes the ring of polynomials in s with coefficients from ¥ q . The group of all 
full rank n x n matrices with entries from ¥ q is denoted by GL(n,¥ q ). The cardinality of 
any set K is given by \K\. The set of positive integers is denoted by Z + . For some integer 
i, we denote the vector in ¥ q , with 1 in the i th position and in the remaining positions, by 
e™. For any matrix M, we denote the submatrix of M, where the row indices run from a to 
b and the column indices run from c to d, by M(a : b, c : d). We denote the a-th row of a 
matrix by M(a, :). A matrix with row vectors v±, V2, ■ ■ ■ , v n is represented as [v±; V2] ■ ■ ■ ; v n ]. 
The column span of a matrix M is denoted by colspan(M). 

2 Multisequences 

We define a sequence S in ¥ q as map from Z to ¥ q . A sequence S = {S(k)}kez, in a finite 
field ¥ q is called periodic if there exists an integer r such that S(k + r) = S(k) for all k. The 
smallest such nonnegative r is called the period of the sequence. There are linear relations 
amongst the elements of a periodic sequence. One obvious example of such a relation being 
S(k + r) = S(k). A general form of such a relation is 

S(k + n) = a n ^iS{k + n — 1) + a n -2S{k + n — 2) + • • • + aoS(k) Vk where a, G ¥ q . (1) 

These are called Linear Recurring Relations (LRRs). The integer n in equation ([T]) is 
called the order of the LRR. Given an LRR we can uniquely associate a monic polynomial 
with it. For example the polynomial associated with the LRR in equation ([1]) is p(s) = 
s n — a n _is n_1 — a n _2<s n ~ 2 — • • • — Q>o- Since we are dealing with periodic sequences, without 
loss of generality, we can assume that 7^ [HJ Theorem 6.11]. 

It is easy to check that all polynomials associated with LRRs of a given sequence S, form 
an ideal 3s in the polynomial ring ¥ q [s}. Since ¥ q [s] is a principal ideal domain, every ideal 
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has a unique monic generating polynomial. The generating polynomial of 3s is called the 
minimal polynomial of the sequence S. The degree of the minimal polynomial is called the 
linear complexity of the sequence. 

Given an LRR of degree n, there are many sequences that satisfy this relation. In fact, 
the collection of all sequences that satisfy this relation form a vector space over ¥ q . The 
maximum possible period of sequences in this vector space is equal to the order of the 
polynomial associated with the LRR. In particular, if the polynomial associated with the 
LRR is a primitive polynomial of degree n, then every nonzero sequence in the corresponding 
vector space has a period equal to q n — 1 ( [TTJ Theorem 6.33]). 

Consider a sequence of linear complexity n. Given n consecutive elements of the sequence, 
every subsequent element can be generated using the LRR corresponding to the minimal 
polynomial. The vector consisting of n consecutive elements of the sequence is called the 
state vector of the sequence. We denote the z-th state vector of the sequence by x(i) i.e., 
x(i) = [S(i), S(i + 1), . . . , S(i + n — 1)]. Observe that if the minimal polynomial of the 
sequence is primitive then the sequence has q n — 1 different state vectors i.e., every nonzero 
vector in F™ is a state vector of the sequence. 

Let aS denote the sequence got by shifting the sequence S once to the left i.e., aS(k) = 
S(k + 1). The k-th state vector of aS is denoted by ax(k). Therefore ax(k) = x{k + 1). 
Observe that ax{k) = x{k + 1) = x(k)A, where 
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This matrix is the companion matrix of the polynomial p(s) = s n — (Ylo~ a i sl )- Observe 
that the companion matrix associated to the polynomial is unique. The sequence obtained 
by shifting S I times to the left, is denoted by a e S i.e a l S{k) = S(k + €). 

Similar to sequences, we define a multisequence in F™ as a map from Z to F™. A 
multisequence W = {W{k)}k& is called periodic if exists a finite integer r such that W{k + 
r) = W(k), for all k. As in the case of scalar sequences, there exist linear recurring relations 
between the elements of the multisequence. These relations are of the form 

W(k + n) = a n _iW(k + n - 1) + a n - 2 W(k + n - 2) H h a W(k) Vk where a; G ¥ q (2) 
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Analogous to scalar sequences, one can define a polynomial p(s) = s n — a n _is n ~ l — a n _2 sn ~ 2 — 
■ ■ ■ — a , which can be associated with every LRR. Again, the polynomials associated to all 
LRRs of a given periodic multisequence, form an ideal in the principal ideal domain F q [s] and 
the monic generator of this ideal is called the minimal polynomial of the multisequence. The 
degree of the minimal polynomial is defined as the linear complexity of the multisequence. 

The i th component of each vector in W gives a sequence of scalars in ¥ q . We call this 
sequence the i th component sequence of W, denoted by W{. Clearly, the minimal polynomial 
of the multisequence is the least common multiple of the minimal polynomials of the com- 
ponent sequences. Therefore, the minimal polynomials of each of the component sequences 
divide the minimal polynomial of the multisequence. Hence, if the minimal polynomial of the 
multisequence is an irreducible polynomial p(s), each of the nonzero component sequences 
also have p(s) as their minimal polynomial. 

Note that a multisequence, with linear complexity n is completely determined by the first 
n terms (vectors). The state of a multisequence can therefore be thought of as n consecutive 
elements (vectors) of the multisequence. Each state is thus an m x n matrix. Thus, we have 
a sequence of matrix states associated with every multisequence. We denote the k-th matrix 
state of the multisequence W by M w (k), i.e. M w (k) = [W(k),W(k + 1), . . . , W{n + k— 1)]. 
The i th row, denoted by Xi(k), of M\y{k) is the k— th state vector of the component sequence 
Wi. In the following theorem we prove that for a periodic multisequence W, all matrix states 
have the same column span. 

Lemma 2.1. For a periodic multisequence, the column span of the matrix states is an in- 
variant. 

Proof. Consider a periodic multisequence W. It is enough to show that colspan(Mw(k)) = 
colspan(Mw{k + 1)), for any given integer k. Let the minimal polynomial of the multi- 
sequence be p(s) = s n — a n _i.s n ~ l — a n _ 2 s™~ 2 — ■ ■ • — a . Since W(k + n) = a W(k) + 
a 1 W{k + 1) + . . . + a n _{W(k + n - 1), therefore W(k + n) e colspan(M w (k)). Thus, 

colspan[M w {k + 1)) C colspan(M w (k)) 

Since a ^ 0, W(k) = ±(W(k + n) - ai W(k + 1) - a 2 W(k + 2) - ... - a n ^W(k + n - 1)), 
i.e., W(k) G colspan(M w (k + 1)). Hence 

colspan(Mw(k)) C colspan(Mw(k + 1)) 
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Therefore colspan(M w (k + 1)) = colspan(M w {k)). Hence proved. 

□ 

We define the dimension of a multisequence as follows 

Definition 2.2. The dimension of a multisequence W is defined as the rank of its matrix 
states. 

As in the case of scalar sequences, any nonzero multisequence with a primitive minimal 
polynomial p(s) of degree n, has a period of q n — 1. In this paper, we henceforth assume 
that the multisequences considered have primitive minimal polynomials. 

The first problem we address is the following: 

Problem 1. Given a positive integer £ and a primitive polynomial p(s) of degree n, how 
many multisequences of dimension £ exist in F™ with p(s) as its minimal polynomial. ( Clearly 
< I < min(m, n) ) 

Two multisequences are considered the same if they are shifted versions of 
one another, i.e., the multisequence W is the same as its shifted version a r W for any r £ Z. 

We denote the collection of £ dimensional subspaces of F™ by G(£, m, ¥ q ). The cardinality 
of G(£, m, W q ) is given by 

Lemma 2.3. Given a primitive polynomial p(s) of degree n, the number of multisequences 
in F™, with minimal polynomial p(s) , having dimension £ is |G(^, m, Fq)| x (q n — q)(q n — 
q 2 )...(q n -q e ~ 1 ). 



Proof. Given a multisequence W of dimension £, by Lemma I2.1[ the column space of the 
matrix states M\y{k) is a unique £ dimensional subspace of F™. Observe that there are 
\G(£, m, ¥ q )\ subspaces of F™ that have dimension £. Consider any one such £— dimensional 
space V. Fix a basis for V, say V\,V2, ■ ■ ■ , v#, where i>j G F™. Let T be the matrix T = 
[vi,V2, ■ ■ ■ ,vg\. Any M £ ]F™ xn w h ose column span is V can then be written as M = TB, 
where B £ ¥ e q xn . The number of such matrices B is equal to (q n — l)(q n — q) . . . (q n — 
q e_1 ) (choosing £ independent vectors in F"). As the polynomial p(s) is primitive, each 
multisequence has q n — 1 distinct matrix states. Thus, the number of multisequences with 
column span V is equal to — — — „_ x — - — - = (q n — q)(q n — q ) . . . (q n — q )■ Therefore, 
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given a primitive polynomial p(s) of degree n, the number of multisequences in F™ with 
minimal polynomial p(s) having dimension £ is \G(£,m,¥ q )\ x (g™ — q){q n — q 2 ) . . . (q n — 
q^ 1 ). □ 

If a multisequence in F™ has dimension m, its component sequences are linearly inde- 
pendent. We can therefore give the following corollary to Lemma [2.31 

Corollary 2.4. Given a primitive minimal polynomial p(s) of degree n, the number of mul- 
tisequences in F™, with minimal polynomial p(s), having linearly independent component 
sequences is (q n — q)(q n — q 2 ) . . . (q n — g m_1 ). 



3 Extensions of Multisequences 

We now look to extend an m-dimensional multisequence W in F™ to an r-dimensional mul- 
tisequence in ¥ r where r > m. Further, we impose a condition that the minimal polynomial 
of the new multisequence is the same as the minimal polynomial of W . An obvious way of 
keeping the minimal polynomial unchanged is by appending to W its component sequences 
or their linear combinations. Thus Wj = Y^iLi a i^i f° r J > m i where € ¥ q . The ex- 
tended multisequence however continues to have dimension m. On other hand, appending 
W with shifted versions of the component sequences may perhaps increase the dimension of 
the multisequence. 

Let R = (fi, r m ) e , with ^2 r& = r. We define the /^-extension of the multisequence 
W in ¥™ as the multisequence Wr in ¥ r q , whose component sequences are obtained from the 
component sequences of W in the following order : Wi, aW\, ...a ri ~ x W\, W2, &W2, ...cr r2_1 W / 2 
, ...Wi, aWi, ...a Ti ~ l Wi, ...W m , aW m , ...a r ' m ~ l W m . Clearly, the minimal polynomial of the mul- 
tisequences Wr and W are the same. We can therefore ask the following question. 

Problem 2. Given R = (7*1, r2, . . . , r m ) 6 Z™, with rj. = r, how many multisequences W 
of rank m in ¥ r q a give R-extended multisequences in ¥ q having dimension r ? 

The solution to this problem is given by the following theorem. 

Theorem 3.1. Let R = (ri,T2, ■ ■ ■ , r m ) e Z™ such that r = ^2 r % an d ^ p{ s ) be a primitive 
polynomial of degree n. The number of multisequences in ¥ 1 q n with minimal polynomial p(s) 
whose R— extensions have dimension r is equal to (q n — q r ~ m+l )(q n — qr-m+2^ _ _ _ _ q r ~ 1 ^ 
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In the remainder of this section we give a constructive proof to this theorem. Starting 
with a multisequence in F™ with dimension m, we recursively generate a series of r — m mul- 
tisequences in F™ culminating in a desired multisequence whose R— extension has dimension 
r. We first prove a few preparatory results which when put together gives us the constructive 
proof. 

For any G = (g\, g 2 , ■ ■ ■ , g m ) G Z™, let G max = maxigi. Let $ define the following map 
from Z™ to Z™. 

$(01.02, • • • ,0m) = (01,02, • • -,90-1,90 - l,0c+l, ■ • ■ ,0m) 

where c is the smallest integer such that g c = G maa ; 

Note that the repeated action of $ on any element of Z+ eventually gives 1 = (1,1, ...,1). 
Thus, given R = (n, r 2 , . . . , r m ) G Z™, $ defines a unique path from i? to 1. We call this 
path the 'R— road'. 

Example 3.2. The R- road for R = (3,2,5,4,1) is (3,2,4,4,1) (3,2,3,4,1) (3,2,3,3,1) 
(2,2,3,3,1) (2,2,2,3,1) (2,2,2,2,1) (1,2,2,2,1) (1,1,2,2,1) (1,1,1,2,1) (1,1,1,1,1). 

Clearly given any point G = (<7i, #2, • • • ,0m) on an R— road, for any other point Q = 
(qi, q 2 , . . . , q m ) lying on the path from R to G, > g, L Vz. Besides, the map $ ensures the 
following: 

- If i < j, gi > gj if and only if gi > rj 

We now look to retrace the i?-road from 1 to R. As a first step, we prove the following 
lemma. 

Lemma 3.3. For every point G = (g±, g 2 , ■ ■ ■ , g m ) ^ R on the R—road, there exists a 
coordinate g c which satisfies at least one of the following conditions: 

1- 9c = G max - 1 and g c < r c . 

2- g c = Gmax and g c < r c . 

Proof. For every point G — (g\, g 2 , ■ ■ ■ , g m ) 7^ R on the R— road, there exists a unique point 
G* on the R— road such that $(G*) = G. Now, G* = (g\, g 2 , . . . , g c -i, g c + 1, 0c+i, • • • , 9m), 
where g c + 1 > Vi 7^ c. Also, since G* is on the path from R to 1, g c + 1 < r c . Therefore, 
g c < r c . If g c + 1 > Vz 7^ c then ^ c = G maa; . If instead, there exists an i such that 
<7 C + 1 = g^ then g c = G max — 1. Hence proved. □ 
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We therefore have the following definition. 

Definition 3.4. Consider an R = (r l7 r 2 , . . . , r m ) G Z™. For every point G = (gi, g 2 , . . . , g m ) 7^ 
R, on the R—road the active coordinate is defined as follows: 

1. If there exists a coordinate g c such that g c = G max — 1 and g c < r c , then the active 
coordinate is the coordinate corresponding to the largest such c. 

2. In the event of there being no g c that satisfies point 1, the active coordinate is the 
coordinate corresponding to the largest c such that g c = G max o,nd g c <r c . 

It can be easily seen that one can traverse the R— road backwards from 1 to R by 
repeatedly incrementing the active coordinate at every point. This is demonstrated in the 
following example: 

Example 3.5. Let R = (3,2,5,4,1). Starting from 1 the R—road is traversed backwards 
as follows: (At every point the active coordinate is underlined) (1,1, 1,1, 1) (1,1,1,2,1) 



(1,1,2,2,1) (1,2,2,2,1) (2,2,2,2,1) (2,2,2,3,1) (2,2,3,3,1) (3, 2, 3,3, 1) (3, 2, 3, 4, 1) (3,2,4,4,1) 



Detecting the active coordinate of any point G involves the following steps: 

• Find G max . 

• Find the largest i such that the i— th coordinate has value G max — 1 and is less than 

n. 

• In the event of there being no i satisfying the preceding condition, find the largest j 
such that the j— th coordinate has value G max and is less than r,-. 

Notice that each of the above steps can be implemented in 0(m) operations. 

For generating multisequences with /^-extensions having maximum dimension, we travel 
backwards along the i?-road from 1 to R. During this backward traversal, at every point 
G on the i?-road, we recursively generate a multisequence whose G-extension has maximum 
dimension. 

We now make the following observation: Given a matrix A G ¥ e q xe in the companion form 
and a vector x — (b±, 62, ■ ■ ■ , be) G F^, for k < £, xA k has the following form 



(3,2,5,4,1). 




k entries 
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where the *s are elements in ¥ q , whose values depend on the matrix A. Therefore, the matrix 
[x; xA; . . . ; xA k ~ x ] has the following structure. 

bi b 2 ... bi^k+i h_ k+2 . . . bi 
b 2 b 3 ... be-k+2 b e _ k+3 ... b e * 

b k b k+ i ... b e * ... * * 

For any G £ Z"\ let N(G, k) denote the number of multisequences in F™ with a given 
primitive minimal polynomial of degree k, whose G— extensions have maximum dimension. 

Theorem 3.6. Let R = (ri, r 2 , . . . , r m ), and let G = (gi,g 2 , ■ ■ ■ , g m ) and be consecutive 
points on the R—road. Then, 

N{G,k) =q m - 1 N(<t>(G),k- 1) (4) 
where k is any integer greater than g = YliLi 9m- 

Proof. Let c be the smallest integer such that g c = G max . Therefore <£>(GQ = (gi, g 2 , ■ ■ ■ , g c -i, g c — 
l,g c+ i, . . . ,g m )- Let W be a multisequence in whose minimal polynomial p k _i(s) is a 
primitive polynomial of degree k — 1. Further assume that the extension of W has 

dimension g — 1. Each matrix state of W is therefore a matrix in ]p™ x ( fc ^ with full row 
rank. As p k _i(s) is a primitive polynomial of degree k — 1, there is a matrix state M of VF, 
whose c— th row is e^l} = (0, 0, 0, 1). For i ^ c, let Xj = [6^, b i2 , . . . , be the z— th 

row of this M. Therefore, M = [xi; x 2 ; . . . ; x c _i; e^lj; x c+1 ; . . . ; x m ]. Now expand M to a 
matrix M* £ F" ixfc as follows: 

1. For every i ^ c, append the z— th row of M with any element of ¥ q . Therefore, the 
i— th row of M* is x* = (x i; d 4 ) £ F^, for some £ ¥ q . 

2. Let the c-th row of M* be i.e., (0, 0, ... , 0, 1). 

Let p k (s) be a primitive polynomial of degree k. Using M* as a matrix state, one can 
generate a multisequence VT* whose minimal polynomial is p k (s). We claim that W* has a 
G— extension with dimension g. 
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As M is a matrix state of W, the following matrix M$(c) is a matrix state of the <3>(G) — 
extension of W: 



M, 



*(G) 



_i ; ■ ■ • ; %c-i] %c-iAk-i] • • • ; x c -iA l }_ 1 ; 

„k-l.k-lA . . „k-l Ag c -2 A . /iSfc+i)- 1 . .„.„/! 

1> e fc— 1 fc— 1' ' ' ' ) e fc— 1 fc— 1 ' -^c+li ^c+l^fc— 1] • • • i ^c+l^fc—i j • • ■ j ^m; J/ m-< 1 fc— 1; • - • j 

^m-rt-k-l J 

where ^4fc_i is the companion matrix of the polynomial pk-i(s) 



The c— th block of rows of M$(g) has the following structure: 

















•• 
•• 



1 

1 * 



£ fl?(Sc-l)x(fc-l) 



For 1 < z 7^ c < m, let re, = {pn , bi2, ■ ■ ■ , b^k-i))- The corresponding z— th block of rows 
of M$(g) has the following structure: 



&i2 



&i2 
bi3 



bi(k-g c ) 
'\k-g c +l) 



bi ai h(g,+l) 



i(k-g c +gi-l) 



bi(k-g c +l) 
h( k -g c+ 2) 



i{k-g c +gi) 



H{k-gi+l) 
\(k- gi +2) 



bi(k-i) 



i(fc-2) 



The *s shown in the blocks above represent entries from ¥ q which depend on the matrix 
Ak-i- Since g c > </j Vz, the *s appear only in the last g c — 1 columns of M$(g) (shown as the 
trailing sub matrix after the vertical line). As $((7)— extension of has rank g — 1, therefore 
M$(G) has rank g—1. 

Similarly, corresponding to the matrix state M* of W*, we have the following matrix 
state of the G— extension of W* . 



M* r * * * * — 1 * * /I * A 92 — 1 * * A * A 

G ~ L X l'i x l^ki ■ ■ ■ ] x l^k ' x 2i x 2 A k] ■■■'■> x 2 A k ! • • • ) X c-1> x c -\^k\ x c -l /i: 

k k A k Age — 1. * * /I * /l3c+l — 1. * * yt 



. ™* /tSm — I"! 
, . . . , ^ m ^L k J 



where is the companion matrix of the polynomial Pk{s). 
For z 7^ c, the z— th block of M G is [a;*; . . . ; (recall that x* = (x i; gQ), where 
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Ak is the companion matrix of the polynomial Pk{s). This block has the following structure 

bil h2 • • ■ bi(k-g c ) 

bi2 h3 ' ' ' fri(fc-g c +l) 

bigi bifa+l) ' ' ' h{k-g c +gi-l) 

The c— th block of M G is [e|; e\A k \ . . . ; e^A^ 0- ]. This block has the following structure. 

••• 
•■■ 

■•• 

Let Mq be the submatrix of M G got by removing its last column and the first row of its 
c— th block. Observe that rank(Mc) = rank(M G ) — 1. By the structure of the c— th block of 
Mq one can clearly see that this submatrix Mq can be modified to M$(g) using elementary 
row operations. Hence this submatrix Mq has rank g—1. This implies that M G has rank g. 
Therefore, W* does have a G— extension with dimension g. 

Note that each of the djS can be chosen in q ways. Each such choice yields a different ma- 
trix M* and hence a different multisequence W*. As a result for every multisequence W 
with minimal polynomial Pk~i{s), the above process gives us q m ~ l multisequences 
W* with minimal polynomial Pk{s). Therefore, 

N(G,k) > g m " 1 iV(<l>(G),A;- 1) (5) 

Conversely, consider a multisequence U* in F™ with primitive minimal polynomial Pk{s) 
whose G— extension has rank g. Consider its matrix state Ml e ]F™ xk whose c— th row is e\. 
Now Ml can be reduced to a matrix M\ G ^ as follows: 

1. For i^c remove the last entry of the i— th row. 

2. Let the c-th row of M l be e k k z\ 

Let Mi generate a multisequence U having primitive minimal polynomial pk-i(s). Using 
similar arguments as those used earlier in the proof, one can prove that the extension 
of U has dimension g — 1. Note that the matrix M\ is independent of the last entries of 



bi{k- 9i +\) bi( k - gt+ 2) 

bi(k- 9i +2) &i(fe- fli +3) 



h(k-i) di 
di * 



bi(k-i) 



•• 
•• 



1 

1 * 



c 9 
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the rows of M*. Hence, there are q m ~ x matrices (including M*), with c— th row e|, which 
have the same first k — 1 columns as Mi. By the above process each one of these matrices 
gives the same matrix Mi (and hence the same multisequence U). Besides if we start with a 
matrix with c— th row e 1 ^ which differs from Mi in any entry corresponding to the first k — 1 
columns, it results in a different Mi (and hence a different multisequence U). Therefore, 

JV($(G),*-1) > (6) 
q m - l N(<S>(G),k- 1) > N(G,k) 

Thus, from equations fl5J) and (jB]) we can conclude that 

iV(G, fc) = g m - 1 A^($(G),A;-l) 



□ 



Using this result, Theorem 13 . 1 1 can by proved in the following manner: 



Proof of Theorem \3.1[ For each j, such that n — r + m < j < n, let Pj(s) be a given 
primitive polynomial of degree j. For every point G = (<7i,<72, • • • ,<7m) on the i?— road, 
let g = XX=i As we have seen in the proof of Theorem 13. 6} starting from a multise- 
quence in F™ with dimension m (i.e., its 1— extension has maximum dimension), having 
minimal polynomial p n - r+m (s) , we can recursively generate multisequences in F™, with min- 
imal polynomial p n - r+g (s), whose G— extensions have maximum dimension, for every G on 
the R— road. 

By Theorem 13.61 for any two consecutive points, 3>(G) and G = (gi,g2, ■ ■ ■ ,g m ) m f ne 
path from 1 = (1,1,..., 1) to R, N(G,n-r + g) = q m ~ 1 N(§(G),n - r + g - 1) where 
g = Y^JiLi g-i- The path from 1 to R has r — m such steps. Therefore, 

N(R,n) = (q m - 1 ) r ~ m N(l,n-r + m) (7) 

However, N(l,n — r + m) is the number of multisequences in F™ of dimension m, with a 
given primitive minimal polynomial p n _ r+m (s) of degree n — r + m. Therefore, by Corollary 
E21 N(l, n-r + m) = (q n ~ r + m - q)(q n ~ r + m - q 2 ) . . . ( q n ~ r + m - q m ~ x ). Hence, 

N{R, U) = ^m-ly-m^n-r+rn _ q ^ q n-r+rn _ g 2y _ _ ( q n-r+m _ q m-lj 

= (q n - q r - m+1 ){q n - q r ~ m + 2 ) . . . ( q n - q^ 1 ) 
Hence proved. □ 
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Remark 3.7. N(R,n) does not depend on the integers (r 1; r 2 , . . . ,r m ) but just their sum. 

One can therefore ask the following question. 

Problem 3. Given any r > m, how many multisequences in ¥ r having dimension r are 
R— extensions of multisequences in F™ for some R = (r 1 ,r 2 , ■ . ■ , r m ) G Z™ where Yl r i = r ■ 

This problem is answered in the following lemma. 

Lemma 3.8. The number of multisequences in which are R— extensions of multisequences 
in F™ is given by, 

N r = ( r ~ 1 W - q r - m+1 )(q n - q r - m+2 ) . . . (q n - (8) 
\r — mj 

Proof. For any r G Z + , define the following subset £R r of Z+. 

m 

SR r :={(r 1 ,r 2 ,..,r m )eZ; | ^r, = r} 

i=l 

Therefore, 

N r = \%.\ x (g n - q r ~ m+1 )(q n - q r ~ m + 2 ) ...(q n - q 1 - 1 ) 

Corresponding to each element of 9^, say (ri, r 2 , . . . , r m ), we can define a monomial, 
Therefore, calculating |$H r | is equivalent to finding the number of monomials of degree r where 
each variable is raised to a nonzero index. However, every such monomial can be written as 
X\x 2 . . . x m T(xi, x 2 , ■ ■ ■ , x m ), where T(xi, x 2 , ■ ■ ■ , x m ) is a monomial of degree r — m. Con- 
sequently, the cardinality of 9l r is equal to the number of monomials of degree r — m. This 
number is equal to f(r-m)+m-i\ = / r -i\ Ag result 

^ V r— m / \r—m) ' 

iV r = ( r ~ 1 ] (q n - q r ~ m+l )(q n - q r ~ m + 2 ) ...(q n - q^ 1 ) 
\r — m J 

□ 

Given R = (r 1 , r 2 , . . . , r m ) G Z™, let r = YnU r i- Let {Pj ( s ) }j=n~r+m be a series of 
primitive polynomials where the index j denotes the degree of the respective polynomial. 
Let AjS be their corresponding companion matrices. Let $(G) and G be consecutive points 
on the R— road. Further, let c be the position of the active coordinate of $(G). Consider 
a multisequence U in F™ with a minimal polynomial p n - r+g -i(s), whose extension 
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has maximum dimension. Note that U can be uniquely determined by any of its matrix 
states. Let My be its matrix state whose c— th row is e^ZlXg-l- The proof of Theorem 13.61 
gives us a procedure to go from My, to the matrix state M v of a multisequence U* with a 
primitive minimal polynomial p n _ r+g (s) of degree n—r+g, whose G— extension has maximum 
dimension. Thus, using this procedure one can generate a sequence of matrices {Mj }Ji n _, r+m 
starting with a matrix M„_ r+m G jp™ x ( n_r+m ) having full row rank, and culminating in a 
matrix M n G ¥ q nx ^ n \ Each matrix Mj in the above sequence uniquely corresponds to a 
point G on the R— road and can be seen as a matrix state of a multisequence with minimal 
polynomial Pj(s) whose corresponding G— extension has maximum dimension. The following 
is an algorithm to generate this sequence: 

Algorithm 3.9. (The variable M is used to store the respective matrix state at every step 
of the algorithm. The current point in the path from 1 to R is stored in the variable G = 
{9ii92i • • • j 9m)- The variable c stores the position of the active coordinate of G. The variable 
g stores the summation of the values of the coordinates of G. ) 
Initialization: 

• Initialize G to 1. 

• Initialize the value of g to m. 

• Initialize M to any matrix in jp™ x ( n - r + m ) /j as j u \\ row rank. 
Main Loop: 

• While g < r do the following 

— Find the position of the active coordinate of G and store it in c. 

D n-r+g 
J n—r+g • 



Find a polynomial f(s) such that M(c, :)f(A n _ r+g ) = e™ 



M = Mf(A n _ r+g ). (This gives us the matrix state whose c—th row is e™_^^. 

For all i ^ c, append the i—th row of M with any di G ¥ q to get the row vector 
(M(i, :),di). Change the c—th row of M to e™Z%+gXx- 

Increment g and g c by 1. 



The most complex part of the above algorithm is to find the polynomial f(s). This can 
be done using the following subroutine. 
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Subroutine 3.10. • Construct the matrix j$ = [M(c, :); M(c, :)A n _ r+g ; . . . ; M(c, :)A™_ r r + g x ] 

• Solve the set of linear equations 

a Jt = e n n -_ r r + + 9 g f° r a e ¥ V +9 - ( 9 ) 

• If a = (ciq, ax, ■ ■ ■ , a„_ r + g _i) zs the solution to the above set of equations, aoM(c,: 
) + axM(c, :)A n ^ r+g + ■■■ + fl n _ r+s _iM(c,:)^:!^" 1 = e^+g- Therefore f(s) = 
a + axs H + a n _ r . +fl _is n ~ r+9-1 . 

Let Ci and C2 be the active coordinates of 1 and respectively. Algorithm 13.91 can be 
thought of as a map from the space of matrices in jp , ^ ix ( n_r '+ m ) which have full row rank and 
whose c\— th rows are e"l£+m) to the space of matrices in F™ xn which have full row rank 
and whose C2— th rows are e™. There are precisely (g n ~ r+m — q)( y q n ~ r+m — q 2 ) . . . (q n ~ r + m _ 
g m_1 ) matrices in flr™ x ( n r+m ) whose C\— th row is e™Ir+m- During each iteration of the 
while loop one can chose rfjS in g m_1 ways. Therefore, corresponding to each choice of 
matrix M n ,_ r+m e F ™ x ("- r + m ) there are g (m-i)(r-m) possi bl e candidates for M n e F™ xn . 
No two distinct choices for the matrix M n ,_ r+m can give the same M n . Therefore, we have 
(q n — q r ~ m+1 )(q n — g r_m + 2 ) . . . (g™ — q r ~ l ) possible matrices which can occur as an output 
to Algorithm 13.91 The number of full row rank matrices in F™ xn , whose Ci~ th row is e™, 
is however (q n — q){q n — q 2 ) . . . (q n — q m ~ l ). Out of these matrices, precisely those matrices 
that occur as matrix states of multisequences whose R— extensions have full rank are the 
ones that can be obtained from the above algorithm. 

We now determine the computational complexity of Algorithm 13.91 We begin by evalu- 
ating the computational complexity of calculating Mf(A n _ r+g ): 

• For any i, M(i, :)A J n _ r+g = (M(i, :)A ] ~^ r+g )A n _ r+g . Therefore, knowing (M(i, :)A^} r+g ), 
M(i, :)A J n _ r+g can be calculated in 0{n 2 ) steps. Consequently, the matrix ^# can be 
generated in 0(n 3 ) operations. 

• Using LU— decomposition, solving the set of linear equations © takes 0(n 3 ) opera- 
tions. 

• For any i, such that (1 < z < n — r + g), the i— th row of Mf(A n _ r+g ) is given by 
aoM(i, :) + axM(i, :)A n _ r+g + ■ ■ ■ + o n _ r+9 _iM(i, As we have already seen 
each element in the above summation can be generated in 0(n 2 ) steps. Each row of 
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Mf(A n _ r+g ) can thus be calculated in 0(n 3 ) steps. Therefore, Mf(A n _ r+g ) can be 
calculated in 0(n 4 ) operations. 

As we have already seen, the active coordinate of G can be found in 0(m) steps. Also, 
appending m — 1 rows of M takes 0(m) operations. Further, incrementing g and g c has 
complexity 0(1). Therefore each iteration of the while loop takes 0(m + n 4 ) time. Since the 
number of steps from 1 to R is (r— to), the while loop runs for a maximum of r — m iterations. 
As a result, the computational complexity of algorithm 13.91 is 0((r — m)(m + n 4 )). 
Therefore for a fixed m and r this computational complexity is 0(n 4 ). This is therefore a 
polynomial time algorithm. 

We now proceed to see an application of the above developed theory. 

4 Word Based Linear Feedback Shift Registers 

The theory developed in the preceding sections finds an application in word based Linear 
Feedback Shift Register (LFSR) design. We begin our discussion by giving a brief introduc- 
tion to Linear Feedback Shift Registers (LFSR)s. 

LFSRs are electronic circuits that implement LRRs. These are widely used in the 
field of pseudo-random number generation and coding theory. LFSRs consist of delay ele- 
ments, feedback elements and adders. For example, the LFSR corresponding to the LRR 




Output I I I I I I I I I I I 

Figure 1: Linear Feedback Shift Register 

S(k + n) = a n ^iS(k + n — 1) + a n -2S(k + n — 2) + • • ■ + aoS(k) is as shown in Figure 1. 
LFSRs with primitive characteristic polynomials are of particular interest since they gener- 
ate sequences with desirable randomness properties like 2-level autocorrelation property and 
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span-n property (all nonzero subsequences of length n occur once in every period) [T]. An 
LFSR can be seen as a state machine where the states are the outputs of the delay blocks. 
Its state transition matrix is the companion matrix of the characteristic polynomial of the 
LRR, (i.e., matrix A in equation (J2])). 

Conventional LFSRs use bitwise operations and hence are incapable of efficiently utilizing 
the parallelism provided by word based processors. In the 1994 conference on fast software 
encryption, a challenge was set forth to design LFSR's which exploit the parallelism offered 
by the word oriented operations of modern processors [T2j . A special case of this scheme was 
implemented by Tsaban and Vishne in their paper [13]. Here, they introduced a family of 
efficient word oriented LFSRs with multiple input multiple output delay blocks. The design 
of Tsaban and Vishne was further generalized in [H] wherein the structure shown in figure 
[2] was proposed to implement the mathematical scheme proposed in |10j . 




Figure 2: Linear Feedback Shift Register with m-Input m-Output Delay Blocks 



Consider the LFSR shown in Figure |2] Let W{k) be the output of the LFSR at the k-th 
time instant. Due to the structure of the LFSR, the following algebraic relation is satisfied 
by the vectors generated by it. 

W(k + b) = B W(k) + B{W(k + 1) + • • • + B b ^W(k + 6-1) (10) 

where e F™ xm . Therefore, for all i, 



W(k + 1) 




W(k) 


W(k + 2) 




W(k + 1) 


W(k + b) 




W(k + b- 1 



(11) 
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where, 
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■mb 
















Bq Bi B 2 



We henceforth call the structure of the matrix A m t> as the m— companion structure. The 
matrix A m b is called the transition matrix of the LFSR and it uniquely characterizes the 
LFSR. The characteristic polynomial of an LFSR is the characteristic polynomial of the 
respective transition matrix. As in the scalar case, LFSRs with primitive characteristic 
polynomials are of special interest. Interestingly, in this design, different combinations of 
feedback matrices can be used to get the same characteristic polynomial. This gives rise to 
the following question: 

Problem 4. Given integers n,m and a primitive polynomial p(s), of degree n = mb, how 
many different LFSR realizations, using m-input m-output delay elements, have p(s) as their 
characteristic polynomial? 

From the discussion above, it is clear that this number is equal to the number of m- 
companion matrices that have the given primitive polynomial as their characteristic polyno- 
mial. Therefore, problem H] can be restated as follows. 

Problem 5. Given integers n,m and a primitive polynomial p(s), of degree n = mb, how 
many m— companion matrices in F™ xn have p(s) as their characteristic polynomial? 

This question was addressed by us in [T3] where the solution was found for the cases 
m = 1, m = 2 and m = n. We will now demonstrate how the theory developed earlier in 
the paper can be used to solve this problem for the general case and in addition give an 
algorithm for generating such configurations. 

Consider an LFSR with m— input m— output delay blocks with primitive characteristic 
polynomial p(s) of degree n = mb. Let A m b be the transition matrix of the LFSR. The 
output of such an LFSR can be seen as a multisequence W = {W(k)}kez in W m . Let 
Wi,W 2 , ■ ■ ■ ,W m be the component sequences of W. Each component sequence of W sat- 
isfies the LRR corresponding to p(s). Therefore, the minimal polynomial of W divides 
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p(s). Since p(s) is a primitive polynomial, the minimal polynomial of W is p(s). The ma- 
trix state of W at time instant k is M(k) = [W(k), W(k + 1), . . . , W(k + n - 1)]. Let 
Xi(k) be the state vector of the component sequence Wi at time instant k. Therefore 
M(k) can also be written as [xi(k); x 2 (k); . . . ; x m (k)]. Let R = (b,b,...,b) and let Wr 

m— times 

be the R— extension of W. Therefore, the matrix state of Wr at time instant k is Mr(A;) = 



[xi(/c); crxi(A;); 



(T 



6-1 



xi(k); x 2 (k); ax 2 (k); 



■ ; a 



b-i 



x 2 {k),...,x m (k)\ax 1 (k)] 



a 



b-l 



x m(.k)\ . 



By permuting the rows of M R (k) we get the matrix Mi (A;) = [xi(k), x 2 (k), . . . , x m (k), axi(k), ax 2 
<jx m (k), . . . , o" 6_1 Xi(/c), <t 6_1 x 2 (/c), . . . , cr b_1 x m (A;)] which can also be written as follows: 



Mi(k) 



W{k) 
W(k + 1) 

W(k + b- 

W(k) 
W(k + 1) 

W(k + b- 



W(k+1) 
W(k + 2) 



W{k + mb-l] 
W(k + mb) 



A 



W(k + b) ■■■ W(k + b + mb-2) 

( W{k) \ ( 
W(k+1) 



rub 



n-l 
rnb 



\W(k + b-l) J 



W{k) \ 
W(k + 1) 



\W(k 



Since Ami, has a primitive characteristic polynomial, for any non zero vector v G F™, the 
vectors v, A m bV, . . . , A m ^v are linearly independent. Therefore, Mi(fc) has rank n. As a 
consequence, Mr(A;) has rank n i.e., Wr has dimension n. We therefore have the following 
lemma: 

Lemma 4.1. Let R = (b,b, . . . ,b) . Given an LFSR with m— input m— output delay blocks 

m— times 

having a primitive characteristic polynomial p(s) of degree n, the R— extension of any non 
zero multisequence generated by it has an dimension n (i.e., maximum dimension). 

Further, since A m i, has primitive characteristic polynomial, repeated action of on 
any nonzero vector d 6 P will generate all nonzero vectors in F" Therefore, by Equation 
pip , for any initial nonzero state of the LFSR all possible nonzero states of the LFSR will 
be covered. As a result, all multisequences generated by the LFSR are just shifted versions 
of each other. In other words, each LFSR with a primitive characteristic polynomial has a 
unique multisequence associated with it. 
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Conversely, consider a multisequence W* = {W*(k)}kez with primitive minimal polyno- 
mial p(s) of degree n, whose R— extension has dimension n. Let A be the companion matrix 
of p(s) and M^(fc) be the matrix state of W* at instant k. One can construct the following 
full rank matrices M*(k) by permuting the rows of the matrix states of the R— extension of 
W*. 



M* w (k) 




M* w {k)A 




_M^(k)A b -\ 





W*(k) W*(k + l) 
W*{k + 1) W*{k + 2) 

W*(k + b-l) W*(k + b) 



W*(k + mb-l) 
W*{k + mb) 

W*{k + b + mb-2) 



M*{k) (12) 



Clearly, for any k e Z, M*(k + 1) = M*{k)A = M*(k - l)A 2 = ■■■ = M*{0)A k+1 . 
Therefore, M*(k + l)M(k)- 1 = M*(0)A k+1 A- k M*(0)~ 1 = M*(0)AM*(0)-\ Thus M*(k + 
1)M(A;) _1 is independent of k and is a constant matrix for a given multisequence W* . Let 
this matrix be denoted by A* mb . Therefore, for all k, 

M*(k + l)=A* mb M*(k). (13) 

Further, given any matrix state M^(k) of W*, A* mb can be constructed as follows: 

A* mh — M*(k)AM*(k)~ 1 . (14) 

where M*{k) is got from Equation (I12p . 

It can be easily verified that the matrix A* mb satisfying equation f[T4|) has an m— canonical 
structure. Let A* mb be as follows: 

" I ... 
7... 

4* - : : : : 

... 7 

R* R* C?* R* 

_±f ^ B 2 ... # 6 _ 

Therefore, for any k, W*{k + b) = B*W*(k) + B{W*{k + 1) + • • ■ + B* b ^W*{k + b-l). Thus 
we have an LRR (and hence an LFSR) generating the the multisequence W. 

Thus, we have demonstrated a one to one correspondence between LFSRs with m— input 
m— output delay blocks having a given primitive minimal polynomial p(s) of degree n = mb 



e F 



mbxmb 
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and multisequences in F™ with minimal polynomial p(s) whose R— extensions have dimension 
n. Therefore, by Theorem 13.11 we have the following: 

Theorem 4.2. Number of LFSRs with m— input m— output delay blocks whose transition 
matrices have a given primitive polynomial p(s) of degree n = mh, as their characteristic 
polynomial is (q n — q n ~ 1 )(q n — q n ~ 2 ) . . . (q n — g n - m+1 ) 

Thus, for R — (b,b, . . . ,b), every multisequence W with a primitive minimal polynomial 

m times 

p(s), such that Wr has dimension mb, is generated by a unique LFSR whose transition 
matrix has characteristic polynomial p(s). Besides, given a matrix state of W one can 
uniquely determine the transition matrix A m b of the LFSR by Equation (fl4l) . 

Therefore, the problem of finding LFSRs generating multisequences with a given primitive 
polynomial reduces to a special case of problem [2] where n = r = mb and R = (b,b, . . .b). 

m times 

Hence algorithm 13.91 can be used to obtain desired LFSR configurations as demonstrated in 
the following example. 



4.1 Example 



We demonstrate Algorithm 13.91 by generating a 3— companion matrix over F 2 with primitive 
characteristic polynomial pe(s) = s 6 + s + 1. We therefore generate a multisequence in F| 
whose (2, 2, 2)— extension has maximum dimension i.e., 6. Note that here R = (2,2,2). 
Consider 



s 5 + s 2 + l 



s 4 + s 



1 
f 1 



which are primitive polynomials over ¥2 of corresponding degrees. Let AiS be the companion 
matrices of the respective Pi(s)s, for 3 < i < 6. We start with a multisequence in F^ with 
minimal polynomial p$(s) and the following matrix state 



1 1 1 
1 1 
1 



We initialize G 
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Iteration 1: G = (1,1,1). Therefore, the active coordinate of G is the 3rd coordinate. (The 3rd 
row of M is already e| and hence M 3 is the desired matrix state). Let us append the 
first and second rows of M 3 with 1 and respectively and change the third row to e\. 
We therefore get the matrix. 

r i 1 1 1 
110 
1 

This is the state matrix of a multisequence W± with characteristic polynomial p±(s). 
Increment the active coordinate of G to get G = (1, 1, 2). It can be verified that Wa is 
a multisequence whose (1,1, 2)— extension has dimension 4. 

Iteration 2: G — (1, 1, 2). Therefore, the active coordinate of G is the 2nd coordinate. The matrix 
state of W4 with second row being e\ is 



Ma 



110 1 
1 
111 



Suppose we append the first and third rows of M 4 with and 1 respectively and change 
the second row to ef . This gives the following matrix: 

110 10 
1 
110 1 

This is the state matrix of a multisequence W5 with characteristic polynomial ps(s). 
Increment the active coordinate of G to get G = (1,2,2). Note that W5 is a multise- 
quence whose (1, 2, 2)— extension has dimension 5. 

Iteration 3: G = (1,2,2). Therefore, the active coordinate of G is the 1st coordinate. The matrix 



state of W5 with first row being ef is 



-Ms 



1 
10 11 
1 1 
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Suppose we append the second and third rows of M 5 with and 1 respectively and 
change the second row to eg. This gives the following matrix: 



Mi 



1 
10 110 
1 1 1 



This is the state matrix of a multisequence W with characteristic polynomial Pq{s). 
Increment the active coordinate of G to get G = (2,2,2). Now W is a multisequence 
whose (2, 2, 2)— extension has maximum dimension, i.e., 6. 

Using the matrix Myy we can construct the following matrix M*: 



1 
10 110 
1 1 1 
1 
10 110 1 
1110 



M* 



M w 
M W A 6 



The 3— companion matrix A 33 can now be obtained as follows: 



A 33 = M*A 6 {M*) 














1 






















1 





1 _ 

















1 




1 





1 


1 










1 


1 





1 





1 




1 


1 


1 


1 





1 



This corresponds to an LFSR whose output multisequence will satisfy the following linear 
recurring recurring relation: 





1 1 




1 




W(k + 2) = 


1 1 


W(k) + 


1 1 


W(k + 1) 




1 1 1 




1 1 
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5 Counting the number of non-singular block Hankel 
matrices 



Hankel matrices are specially structured matrices which frequently appear in the fields of 
signal processing [16], image processing and control theory. In this section we derive a 
formula for the number of non-singular Hankel matrices of a given size, over a given finite 
field W q , by using the theory developed in the preceding sections. 

A Hankel matrix is a matrix which is constant along the anti-diagonals. For example: 



H 



(l\ (12 ... Gtn-1 O n 

O2 O3 • • • a n O n+ i 



On— 1 O n 
On O n +i 



02n-3 02n-2 
02n-2 02n-l 



(15) 



It can be easily seen that the space of n x n Hankel matrices is a 2n — 1 dimensional 
space and each Hankel matrix can be uniquely determined by the corresponding vector 
an = (ai, 02, ■ • ■ , 02 n _i) G Fg' ra_1 . Since the number of Hankel matrices over a finite field ¥ q 
is finite, one may pose the following question. 

Problem 6. Given n, find the number of Hankel matrices in F'" xn that have full rank. 



We solve this problem by proving a bijection between the set of full rank Hankel matrices 

2 



in F" xn 



and the set of multisequences in ¥ q with a given primitive minimal polynomial p(s) 



whose R— extensions have maximum dimension, for R 



[n 



l,n). 



Theorem 5.1. Let p(s) be a primitive polynomial of degree 2n — 1. Let R = (n — l,n). 
Consider a Hankel matrix H G F™ xn corresponding to the vector an = (01, 02, . . . , 02n-i) G 
Fg" -1 . The matrix H has full rank if and only if the matrix M = [e^Zi, o#] is a matrix 
state of a multisequence W with minimal polynomial p(s) , whose R— extension has maximum 
dimension. 



Proof. Consider the multisequence W with primitive polynomial p(s) and matrix state M . 
Therefore the corresponding matrix state Mr of the R— extension of W is as follows: 
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. 














1 





. 











1 


* 





. 








1 


* 


* 


a i 


a 2 


a n _i 




a n+l 


• • • &2n-2 




a 2 


a 3 . 


a n 




0"n+2 


• • • &2n-l 


* 




a n . 


■ «2n-3 


^2n-2 


0>2n-l 


* * 


* 


a n 




■ «2n-2 


^2n-l 


* 


* * 


* 



Clearly, the submatrix Mr(1 :n — 1, n + 1 : 2n — 1) (the top right submatrix) has full rank. 
Therefore, Mr has full rank if and only if the submatrix Mr(u : 2n — 1, 1 : n) (the bottom 
left submatrix) has full rank. However, M#(n : 2n — 1, 1 : n) is the Hankel matrix if. Hence, 
the matrix H has full rank if and only if the matrix Mr has full rank. In other words, the 
matrix H has full rank if and only if the R— extension of the multisequence W has maximum 
dimension. Hence proved. □ 

Every multisequence W in F^, with primitive minimal polynomial p(s) is uniquely char- 
acterized by any of its matrix states. Therefore, the number of full rank Henkel matrices in 
F™ xn is equal to the number of multisequences in Fjj whose R— extensions have maximum 
dimension. Hence, by Theorem 13.11 we have the following theorem 

Theorem 5.2. The number of Hankel matrices in F^ xn having full rank is (q 2n ~ l — q 2n ~ 2 ). 



6 Conclusions 

In this paper we have introduced the concept of matrix states. Using matrix states, we 
have defined the dimension of a multisequence and calculated the number of multisequences 
with a given dimension. The concept of R— extensions has then been introduced. We have 
calculated the number of multisequences whose R— extensions have maximum dimension. 
Further we give an algorithm to generate such multisequences. We have then demonstrated 
an application of the theory developed for R— extensions. For any given m, we have derived 
a formula for the number of LFSR configurations, with m input m output delay blocks, 



25 



that generate multisequences with a given primitive minimal polynomial. Further, we have 
demonstrated the use of the algorithm developed for R— extensions, for the generation of 
such LFSR configurations. Finally, using the theory developed, we have derived a formula 
for the number of Hankel matrices in F™ xn that have full rank. 
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